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Abstract 

We investigate the two-body decay modes pp, 7T7t*(1300) and a a of the /o(1500), 
all leading to the 4-7r decay channel, in a three-state mixing scheme, where the /o(1500) 
is a mixture of the lowest-lying scalar glueball with the nearby isoscalar states of 
the 0 ++ QQ nonet. In the leading order of this scheme, the decay mechanism of the 
/o(1500) proceeds dominantly via its quarkonia components, which can be described in 
the framework of the 3 Po QQ pair creation model. We predict the hierarchy of decay 
branching ratios B with B(pp)^B(Tnr)^B(aa) > B(tttt*), providing a key signature 
for the proposed mixing scheme in this leading order approach. 
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Recently, there has been growing attention, both experimental and theoretical, on the 
possible existence of glueball states. A scalar-isoscalar resonance / 0 (1500) is clearly estab¬ 
lished by Crystal Barrel at LEAR i a in proton-antiproton annihilation, which is regarded 
as a rich gluon production mechanism with favored coupling to glueballs ||. The main 
interest in the / 0 (1500) as a candidate reflecting the intrusion of a glueball state in the 
scalar-isoscalar meson spectrum rests on several phenomenological and theoretical observa¬ 
tions. Lattice QCD, in the quenched approximation Q, predicts the lightest glueball to be a 
scalar, lying in the mass range of 1500-1700 MeV. Although the decay pattern of the /o(1500) 
into two pseudoscalar mesons is compatible with a quarkonium state nn = [uu + dd) in 
the scalar QQ nonet, the observed total width of T ^ 120 MeV is in conflict with naive 
quark model expectations of T > 500 MeV 0. Moreover, it is now believed 10. that the 
/o(1370) state occupies the corresponding position in the scalar nonet, with decay branching 
ratios and total width consistent with a dominant nn structure. As the / 0 (1500) cannot be 
accommodated as a radial excitation of the scalar 3 P 0 nonet [0, it appears as an excess state 
in the QQ systematics. 

In the analysis of the partial decay widths into two pseudoscalar mesons, several schemes 
have been proposed to attribute at least a partial glueball nature to the / 0 (1500) state 

|8|, although quantitative predictions for its glueball content differ sizably. In the 
scenario proposed by Ref. 0 the pure glueball is mixed with the nn and ss states of the 
scalar 3 Po nonet. The mixture leads to three scalar states with nearly equal glueball strength, 
identified with the / 0 (1370), / 0 (1500) and an additional state around 1700 MeV which can 
possibly be the /j=o, 2(1710). It is argued that in the leading order of strong coupling QCD, 
the hadronic two meson decay proceeds dominantly via the QQ- components of the /o(1500), 
if the Fock states of final state isoscalar mesons contain no gluonic component. Neglecting 
an intrinsic violation of flavor symmetry, the model can explain the decay pattern of the 
/o(1500) into two pseudoscalar mesons. In particular, the experimental weak KK decay 
mode is obtained from the destructive interference between the admixed nn and ss modes. 

There are hints that the scalar ss state may have a lower mass than the scalar glueball 
motivating another three-state mixing model with this level ordering [9;]. In this scenario 
the glueball component dominantly resides in the fj( 1710) with the / 0 (1500) remaining as a 
nearly pure quarkonium state, but the destructive interference between the nn and ss flavors 
for the / 0 (1500) still remains. For this reason, both analyses of Refs. || and 0 are very 
similar for the decay of / 0 (1500) into two pseudoscalar mesons 0. 

In the present work, we go beyond the two pseudoscalar meson decay modes and study 
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the predictions of the three-state mixing schemes of Refs. [6;. |9j in the 47 t decay modes. We 
derive the decay pattern of the / 0 (1500) into the two-body decay channels pp, 7T7t*(1300) 
and a a (a is the broad tttt S-wave resonance). They all contribute to the 47 t decay rate, and 
we calculate their relative ratios to that of the 27r decay channel. The two-body decay modes 
leading to the 47r final state have the advantage that, in the decay dynamics, no SU(3 ) flavor 
symmetry breaking and mixing angles for the meson nonets in the decay channel are needed. 
This reduces the uncertainties occuring in the analysis of the two pseudoscalar meson decays. 

The model - As mentioned earlier, the hadronic decay of the admixed /o(1500) is domi¬ 
nated by the QQ components. The coupling of the scalar QQ components to the final 
two-meson states can be calculated in the non-relativistic QQ pair creation 3 P 0 model [|1C)|], 
in which the QQ pairs are created with vacuum quantum numbers. The relevant OZI-allowed 
quark line diagram, describing the flavor flux from the initial QQ component of the /o(1500) 
to the two meson final states BC (=7T7t, pp, aa and 7T7r*(1300)) is shown in the figure. The 
details of all the following formulae can be found in Refs. § 0 m [TT| . The transition 
dynamics is governed by the nonperturbative QQ 3 Pq vertex: 
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p 3 ( 4 ) are the momenta of the created QQ pair and Tf /t (p) = IpI^^p)- The identity op¬ 
erators If and 1 c project onto singlet states in flavor and color space, respectively. The 
dimensionless parameter A corresponds to the strength of the transition, which in turn is 
related to the probability for QQ pair creation out of the hadronic vacuum. Since for the 
considered final states BC no ss pair creation is involved A is flavor independent. With 
harmonic oscillator wave functions for the quark clusters, the scalar QQ-component of the 
/o(1500) is given as: 

$ /o(Pl.P2) = lpS fo =l,l fo =l,J fo =0(Pl,P2)X f FXc 
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Rf 0 is a size parameter, the intrinsic spin \s =1 °f the QQ pair is coupled with relative orbital 
angular momentum l — 1 to the J PC = 0 ++ state, \c denotes the color wave function and 
the effective flavor part Xf gi ven by Xf = £ ( n 'h ), where £ represents the amplitude for the 
nn component in the three-state mixing scheme. In the same way, for the final state mesons 
a and it* (which is the first radial excitation of the 7r with n n * = 1) with decay momentum 
K and intrinsic quark momenta Pi(j), the wave functions are 
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with the Laguerre polynomial = 3/2 — p. Here, a is effectively described by an 

isoscalar 3 P 0 QQ state. The S-wave mesons 7r and p are described by 
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with Sn = J-n = 0, S p = J p = 1. The size parameters for 7r and 7r* have to fulfill the relation 
R n = R n * due to orthogonality of the wave functions. 

The transition amplitude for the process /o(1500) —> BC of the figure is given by 


T = < t2 3) (pi,P3)T^ j (p 4 ,P2)|H 3 ^|T /o (pi,P2) > 


,(42) 


(34) | 


( 6 ) 


_ rpfo^BC rpfo^BC rpfo^BC 

— 1 SS 1 F 1 C 


The color part Tq is 4=, and the spin-spatial part Ts°g^ BC is given with the appropriate wave 
functions of Eqs. (2)-(5) as 
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where the internal quark momenta are integrated over and evaluation of the integral (7) is 
done analytically ||10|| . The various flavor factors Tp°^ BC> are indicated in the table. Using 
relativistic phase space, the partial decay width T f 0 ^BC is written as 
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with decay momentum |K| = K and where = y M 2 + K 2 and M t (i=B,C) are energy 
and mass of the outgoing meson i. Analytical evaluation of Eq. (|8]) leads to 
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with A = 1, i? = 9 for / 0 —> 7T7T and A = 3, B = 1 for f 0 —» pp. For the broad mesons 
/o(1500), 7r*(1300), p(770) and cr(760) we average over the mass spectrum f(m ): 

T fo ^ BC = J dm fo dm B dmc T fo ^ B c(m fo ,m B ,mc)f(rrif 0 )f(m B )f(mc) (12) 
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with a proper threshold cutoff introduced as in Ref. |T2|. Masses Mj and widths E; are taken 
from the Particle Data Group. They are well defined except T^* which is (400 ± 200) MeV. 
We parameterize the mass distribution of the scalar-isoscalar a meson as in Ref. [[p| ( M a = 
760 MeV, To- = 640 MeV). The partial decay widths of all An channels and also the 2n mode 
scale with (£ • A) 2 . The predictions for relative branching ratios are therefore independent of 
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any particular three-state mixing scheme [|6j. j||, specified by the explicit value of £. The only 
parameters left in the model are the size parameters R n , R p , R a and Rf 0 of the meson wave 
functions. Assuming a common oscillator parameter for all mesons, meson decay analyses 
yield R = 2.5 Gel/ -1 || |TlJ . These values will be refered to as set 1. In Ref. @, for the 
two pseudoscalar meson decay modes, differences in size parameters for the QQ states are 
taken into account, with i?( 3 P 0 ) = 2.0 GeR -1 and P( 1 S' 0 ) = 2.5 GeR -1 (set 2). A recent 
full scale analysis of tensor meson decays 0 yields R n = R p = 3.69 GeR 1 , suggesting 
R( 3 Po ) ~ 3.7GeR -1 (set 3). There is no reliable constraint for R a ; since a represents 
a broad tht S-wave resonance, R a > R n is reasonable. We resort to these different sets of 
values to study the sensitivity of the different observables to this input leading to an effective 
error band of the decay model. 

Results - If we adopt R n = R a and T^. = 400 MeV we obtain for the fo —>■ BC decays 
the following ratios: 
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: 0.93 

: 0.33 

with set 1 

1 

: 1.04 

: 0.76 

: 0.56 

with set 2 

1 

: 0.92 

: 0.28 

: 0.19 

with set 3 
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We can see that all decay rates reveal a rather strong dependence on the size parameters, 
ffowever, even when we vary the size parameters in a rather wide range, the 47 t decay modes 
of the / 0 (1500) always follow the decay pattern 


< B(crcr) ^B(tttt) <B(pp) 


(14) 


This hierarchy of the decay modes remains also unaltered when taking into account the 
experimental uncertainty of T^. by ±200MeR. Similarly, a variation of the a size parameter 
in the range R n < R a < 6.0 GeR -1 does not change the order (p~4|). This constitutes a 
stable prediction of the model. Experimental data concerning the 47T decay mode of the 
/o(1500) and its separation into the contributions of the individual two-body channels are 
still incomplete. Recently, Crystal Barrel observed the /o(1500) in its 47r° decay mode in 
the annihilation reaction pp —> 57t° at rest |15|. With their results for / 0 (1500) —> 7T7t, given 
by a single channel (pp —> 37t°) [[I]] and, more consistently, by a coupled channel analysis 
(pp —> 7r°?7?7, 7r°7r°?7 and 37 t°) |lf|, the ratio r is found to be 


r = 


P(/o(1500) -> 4 tt) 
P(/o(1500) -> 2 t r) 


3.4 ±0.8 [15] 
2.1 ±0.6 


(15) 








where the contribution of the pp intermediate state is not included in P(/ 0 (1500) —> 47t). 
However, an alternative multi-channel analysis [|T^], resulting in a partial decay width of 
r= 60 ± 12 MeV and a total width of Y = 132 ± 15 MeV for the / 0 (1500) seems to 
indicate a ratio of r closer to 1. For the three parameter sets and using = 400 MeV, 
B(n* —> an) = 0.6, and R a = R n , our calculation gives 


r = 


Bjjp -> 4tt) 

B(f 0 - 2vr) 


1.13 for set 1 
1.09 for set 2 
0.40 for set 3 


(16) 


to be compared to the experimental number of Eq. ([15|). Again, taking into account the 
experimental uncertainties in and in B{n* —> an) by ±0.2 the results (|16j) are changed 
by less than 10%. A variation of R a from the value of R n to 6 GeV~ l leads to ranges of 
predictions for r of 0.41 — 1.21, 0.32 — 1.20 and 0.26 — 0.45, for the three parameter sets 
respectively. Even when one allows strong variations in the size parameters the theoretical 
predictions for r are lower than the experimental result by at least a factor of two. Although 
in Ref. [Hj, the hnal states a a and nn* contribute with similar intensities to the 47 t° mode, 
precise values for the respective decay branching ratios cannot be given due to interference 
effects. Furthermore, preliminary results |T9|] of an analysis of pp —> 57T including charged 
state combinations of n, seem to indicate that the aa decay is dominating over the pp decay. 

In summary, we have investigated the glueball-quarkonia mixing scheme of Refs. |6|, [| 
in the / 0 (1500) decay into aa, pp, n*n and 7T7r hnal states. The leading order hadronic decay 
mechanism, as described in Ref. |(||, proceeds both by the QQ and the gluonic component of 
the /o(1500); latter process, where the glueball decays into a pair of glueballs, can feed hnal 
states with isoscalar mesons, possessing a sizable gluonic component. Given the modelling 
of a as a 3 P 0 QQ state, where gluonic components are neglected, for the leading order 
decay mechanism we predict a decay pattern of 47r modes, where the pp channel dominates. 
This is in conflict with preliminary experimental results [I9|, where the pp decay mode is 
suppressed. For the 47r/27r decay ratio the result we obtain is insufficient to fully explain 
the experimental result of Eq. (p±|). Another analysis of the / 0 (1500) decays 


seems to 


indicate that the ratio of Eq. (jl5|) is considerably lower, reducing the disagreement with the 
theoretical results. We stress that the decay pattern developed here is common to a simple nn 
state and to the mixed / 0 (1500), where its ss configuration does not couple to the 47r modes. 
If the experimental qualitative ordering of the / 0 (1500) decay rates B{An) > B(2n) and 
B(aa) > B(pp) is confirmed, the decay scheme developed here should be modified by taking 
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into account an additional direct coupling of the gluonic component of the / 0 (1500) state to 
the 47 t decay channels; for example, leading order contributions can enhance the crcr decay 
mode ||, while higher order decay mechanisms can contribute to all 47T decay channels. Thus 
an accurate experimental determination of the An decay modes of the / 0 (1500) provides an 
additional sensitive signal for the need to go beyond a pure QQ consideration of the / o (1500), 
possibly revealing the significant glueball configuration claimed to reside in this state. 
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